An RXTE Archival Search for Coherent X-ray Pulsations in 

LMXB 4U 1820-30 
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ABSTRACT 

As part of a large-scale search for coherent pulsations from LMXBs in the RXTE archive, 
we have completed a detailed series of searches for coherent pulsations of 4U 1820—30 — an 
ultracompact LMXB with a binary period of 11.4 min, located in the globular cluster NGC 6624. 
The short binary period implies any coherent signal would be highly accelerated, so we used 
phase modulation searches, orbital-parameter-fitting coherent searches, and standard acceleration 
searches to give significant sensitivity to millisecond pulsations. We searched, in four energy 
bands and at a range of luminosities, a total of 34 archival RXTE observations, 32 of which 
had on-source integration times longer than 10 ks, and some of which were made consecutively 
which allowed us to combine them. We found no pulsations. Using our phase modulation search 
technique, which we ran on all 34 observations, we have been able to place the first stringent (95% 
confidence) pulsed fraction limits of < 0.8% for all realistic spin frequencies (i.e. < 2 kHz) and 
likely companion masses (0.02 M Q < M c < 0.3 M©). Using our orbital-parameter-fitting coherent 
search, which we ran on only 11 selected observations, we have placed a pulsed fraction limit of 
< 0.3% for spin frequencies < 1.25 kHz and companion masses M c < 0.106 M©. By contrast, all 
five LMXBs known to emit coherent pulsations have intrinsic pulsed fractions in the range 3% 
to 7% when pulsations are observed. Hence, our searches rule out pulsations with significantly 
lower pulsed fractions than those already observed. 



Subject headings: binaries: 
individual(4U 1820-30) 

Introduction 
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One of the great scientific expectations when 
the Rossi X-ray Timing Explorer (RXTE) was 
launched in 1995 was the discovery of coher- 
ent pulsations from low-mass X-ray binaries 
(LMXBs). At present, only five accreting mil- 
lisecond pulsars in LMXBs are known. All five 
are faint transient sources for which the pulsa- 
tions at the spin period of the pulsar were discov- 
ered during outbursts (Wijnands & van der Klis 



1 Department of Physics, McGill University, Montreal, 
QC H3A 2T8 

2 E. O. Hulburt Center for Space Research, Code 7655, 
Naval Research Laboratory, Washington, DC 20375 

3 Department of Mathematics and Statistics , McGill 
University, Montreal, QC H3A 2K6 



1998; Markwardt et al. 2002; Galloway et al. 2002; 
Markwardt, Smith, & Swank 2003a; Markwardt 
& Swank 2003). 

Surprisingly, no pulsations have ever been seen 
in the non-transient X-ray emission from any 
LMXB. One possibility is that MSP coherent pul- 
sations are smeared out by rapidly changing or- 
bital Doppler shifts in short period binaries. Mod- 
ern accelerated signal search techniques have not 
been applied systematically to many LMXBs so it 
is possible that there are millisecond pulsations in 
them that have gone unnoticed. 

There are many things that we can learn from 
searching for more examples of direct pulsations 
from LMXBs. Coherent pulsations give the pre- 
cise rotation rate of the neutron star and test the 
connection between recycled MSPs and LMXBs. 
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They help constrain the models of kHz QPOs and 
burst oscillations. Timing of the coherent pulses 
would also allow measurements of the accretion 
torques, the orbital parameters, and potentially 
the companion masses. In addition, in combina- 
tion with X-ray and optical spectra, pulse tim- 
ing can give information about the compactness 
and the equation of state of the neutron star. Fi- 
nally, setting stringent upper limits on coherent 
pulsations in several sources in a variety of spec- 
tral states will impose constraints on the possible 
mechanism for suppression of the predicted coher- 
ent pulsations mentioned above. For these rea- 
sons, we have started a large-scale archival RXTE 
search for coherent pulsations from LMXBs, the 
first target of which is the ultracompact globular 
cluster LMXB 4U 1820-30. 

This paper is organized as follows. In Sec- 
tion 1.1, we describe the general properties of 
4U 1820—30. In Section 2, we describe the obser- 
vations as well as the procedure we used to search 
for coherent pulsations. In Section 3, we describe 
our search techniques and the methods used to 
investigate the obtained candidate pulsation fre- 
quencies. In Section 4, we describe our results 
and derive our upper limits on possible coherent 
pulsations. Finally, in Sections 5 and 6, we discuss 
the implications of our findings and summarize. 

1.1. Characteristics of 4U 1820 30 

The source 4U 1820-30 is an atoll LMXB in 
globular cluster NGC 6624 (Stella, Pricdhorsky, 
& White 1987). It has an orbital binary period of 
685 s (11.4 min). The orbital period was identi- 
fied to be the period of a nearly sinusoidal X-ray 
luminosity modulation with an amplitude of 2% 
to 3% and no strong energy dependence (Stella 
ct al. 1987; Smale et al. 1987; Sansom et al. 1989; 
Tan et al. 1991; van der Klis et al. 1993; van 
dcr Klis et al. 1993; Chou & Grindlay 2001). A 
16% modulation at a period of 687±2.4 s was 
also seen in the UV band (Anderson et al. 1997). 
This is the shortest known binary orbital period 
in an LMXB. 4U 1820—30 undergoes a regular 
~176 day accretion cycle (Priedhorsky & Terrell 
1984) switching between high and low luminosity 
states, which differ by a factor of ~3 in their lu- 
minosity (Stella et al. 1987). This luminosity vari- 
ation makes 4U 1820—30 an attractive target for 
our search because it allows us to sample a fairly 



large range of accretion rates in case the pulsa- 
tions become visible only in some. In the low 
luminosity state, regular Type I bursts are seen 
±23 days around the minimum luminosity (Chou 
& Grindlay 2001). It is notable that burst oscil- 
lations have not been observed from this source. 
In the low state, 4U 1820—30 has also shown 
an extremely energetic superburst (Strohmayer & 
Brown 2002). This bursting behaviour confirms 
that the observed luminosity variations actually 
correspond to changes in the local accretion rate 
at the surface, and not just a geometrical effect, 
since they affect the bursting properties of the 
source. 4U 1820—30 also exhibits a variety of 
spectral states as can be seen from Figure 1. Sev- 
eral low frequency QPOs (Stella, Pricdhorsky, & 
White 1987; Dotani et al. 1989; Wijnands, van der 
Klis, & Rijkhorst 1999). as well as two peaks of 
kHz QPOs (Smale, Zhang, & White 1997; Zhang 
ct al. 1998; Kaaret et al. 1999) have been observed 
from this source. The separation between the two 
peaks is 275±8 Hz (van der Klis 1999). If the sec- 
ondary star in the system is a white dwarf, its mass 
is estimated to be 0.058 to 0.078 M Q (Rappaport 
et al. 1987). If the secondary is a main sequence 
star, the upper limit on its mass is 0.3 M (Arons 
& King 1993). 

Of all known LMXBs, we searched 4U 1820-30 
first for several reasons. First, many high-time- 
resolution, long data sets are available in the 
archives. Second, the presence of kHz QPOs 
in this source may indicate favorable conditions 
for detecting pulsations. Third, its extremely 
short orbital period allows the use of the phase 
modulation search technique. This technique is 
most sensitive and provides the biggest increase 
in sensitivity over previous methods when the 
observations are longer than two complete bi- 
nary orbits (Ransom, Cordes, & Eikenberry 2003). 
The precisely known and short orbital period also 
makes 4U 1820-30 one of only two LMXBs that 
can be searched using our fully coherent orbital- 
parameter-fitting search technique on a reasonable 
time scale. 

The source 4U 1820—30 has previously been 
searched for coherent pulsations at frequencies un- 
der 50 Hz by Wood et al. (1991) and by Vaughan 
ct al. (1994). Their upper limit on the pulsed 
fraction for sinusoidal pulses was 0.6% (95% con- 
fidence level), and 1.1% (99% confidence level) re- 
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spcctively. Our searches cover a much wider fre- 
quency range, including, for the first time, the kHz 
range. This is particularly important since likely 
spin frequencies are near 275 Hz or 550 Hz, based 
on the kHz QPOs (van der Klis 1999; Chakrabarty 
ct al. 2003). 

2. Observations and Preliminary Data 
Preparation 

We searched 34 archival RXTE PCA observa- 
tions collected between 1996 and 2002 (see Ta- 
ble 1). All the observations were available in event 
modes with time resolutions of ~125 [ms. 32 of the 
observations had a total time on source > 10 ks. 
The longest of these was 29 ks with a total time 
on source of 16 ks. Some of the observations were 
segments of longer ones which allowed us to con- 
catenate and analyse them together. The longest 
of the concatenated observations was 77.5 ks with 
a total time on source of 46.5 ks. 

For each observation, we downloaded the files 
of raw events (one file for each RXTE orbit). The 
duration of an RXTE orbit (including Earth occul- 
tation time) is ~5.7 ks. We filtered the files by se- 
lecting only events where the source elevation was 
> 10°, the pointing offset was < 0.02°, and where 
a minimum number of operational PCUs were on 
(see Table 1). We then transformed the arrival 
times of the events to the solar system barycenter 
using the 'FTOOL' 1 FAXBARY. The barycentering 
was necessary to remove the effects of the space- 
craft and Earth motions. We then split the data 
into four energy bands: Soft (2-5 keV), Medium 
(5-10 keV), Hard (10-20 keV), and Wide (2-20 
keV), in order to attempt to maximize the signal- 
to-noise ratio of the unknown pulsations. This 
processing used custom Python scripts to call the 
standard 'FTOOLs'. At the end of the preliminary 
data preparation we obtained four lists of events 
(photon arrival times) for every RXTE orbit in ev- 
ery observation, each list corresponding to one of 
the four energy bands. 

3. Data Analysis 

After we prepared the RXTE-orbit-long lists of 
events, we concatenated them into observation- 
long lists of events, and also broke them into much 

http://heasarc.gsfc.nasa.gov/ftools/ 



shorter lists of events. We then binned the events 
into time series of different time resolutions (see 
Table 2) and searched them for pulsations using 
three types of searches which we summarize in Sec- 
tions 3.1.1, 3.2.1, and 3.3.1. In Sections 3.1.2, 
3.2.2, and 3.3.2, we explain how we used the three 
searches to obtain candidate pulsation frequencies. 
Finally, in Sections 3.1.3, 3.2.3, and 3.3.3, we ex- 
plain the methods we used to investigate the ob- 
tained lists of candidate pulsation frequencies. All 
the data sets were processed on the 52-node dual 
1.4 GHz Athlon-based Beowulf cluster operated by 
the McGill University Pulsar Group. 

3.1. Acceleration Searches 

3.1.1. The Search Technique 

This is the method traditionally used to com- 
pensate for the effects of orbital motion (e.g., 
Johnston & Kulkarni 1991): a time series is 
stretched or compressed appropriately to ac- 
count for a trial constant frequency derivative 
(i.e. constant acceleration), Fourier transformed, 
and then searched for pulsations. We used a 
Fourier-domain variant of this technique based 
on matched filtering techniques which includes 
Fourier interpolation and harmonic summing of 1, 
2, 4, and 8 harmonics to improve the sensitivity to 
non-sinusoidal pulsations (Ransom, Eikenbcrry, & 
Middleditch 2002). 

Acceleration searches work best for systems 
where the orbital period ~ 10 x the observation 
length. Acceleration searches are not optimal for 
systems like 4U 1820—30 where the orbital period 
is shorter than the observation duration: the ac- 
celeration is not constant over the length of our 
time series. However, we chose not to run searches 
on time series shorter than 600 s (see Table 2) be- 
cause we required a large number of photons. So 
while we did not expect acceleration searches to 
find pulsations, we still used them because they 
are included in our standard search pipeline. Also, 
the segmented acceleration searches that we ran, 
even though they were not short enough for the 
acceleration to be constant, can still be useful if a 
strong pulsed signal appears for a short period of 
time perhaps because the 'weather' on the surface 
briefly cleared and the pulsations appeared. 
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3.1.2. Obtaining Candidate Pulsation Frequen- 
cies 

To run acceleration searches, we first binned the 
lists of events which we obtained in the data prepa- 
ration into time series of different lengths and time 
resolutions (see Table 2). We performed an FFT 
on each of the time series and fed the result into 
a local red noise reduction algorithm which di- 
vides the FFT into continuous intervals the size 
of which increases logarithmically with frequency. 
The power in every bin of every interval is divided 
by the mean power in that interval. This produces 
FFTs that have a reduced signature of red noise. 
We then ran the acceleration searches on the ob- 
tained FFTs, searching for a maximum constant 
acceleration of a max < 5.1 x 10 10 /(.fT 2 obs ) m/s 2 , 
where / is the frequency of the detected pulsations 
in Hz, and T b s is the integration time in seconds. 
For a signal of period 3 ms and an integration 
time of 600 s, this acceleration is 425 m/s 2 . This 
number was chosen because it corresponds to the 
maximum length of a trial Fourier response tem- 
plate used in the search that would fit in the cache 
memory of the available machines. The maximum 
possible line of sight acceleration of 4U 1820—30 
is a max (x) — 25240.6 x x m/s 2 where x is the 
projected semi-major axis of the binary in light- 
seconds. For a companion mass of 0.05 M Q and 
an inclination angle of 60°, x — 13.4 lt-ms, and 
a m ax = 339 m/s 2 . 

Using acceleration searches, we searched for 
signals with an equivalent Gaussian significance 
a > 2 (Groth 1975) and with frequencies between 
0.0375 Hz and 2048 Hz. We estimated a based on 
the probability distribution of sums of powers in 
a normalized power spectrum (Groth 1975; Ran- 
som et al. 2002). In order to correctly determine 
the significance of a candidate, we also took into 
account the number of independent trials in the 
searches (Ransom et al. 2002). Sometimes we in- 
creased the lower limit on the searched frequency 
to 1/8 Hz in order to eliminate the large number of 
low-frequency candidate oscillations that the red 
noise reduction routine did not successfully elimi- 
nate. 

3.1.3. Candidate Follow-up 

When our acceleration searches detect a can- 
didate oscillation, they return values for the fre- 



quency of the signal /, the frequency derivative /, 
and the significance of the detection a. Using the 
returned values for / and /, we folded all acceler- 
ation candidates with a > 3. Our folding software 
had the ability to fine tune the search in the / 
and / directions. We also folded all the candidates 
with a > 2 that appeared in 15 or more observa- 
tions to within 0.1% of their frequency. 0.1% was 
chosen because the frequency modulation caused 
by the the orbit of 4U 1820-30 will not cause / to 
deviate by more than 0.1% (see Section 3.2.1). All 
the folded candidates showed noise-like profiles. 

3.2. Phase Modulation Searches 

3.2.1. The Search Technique 

In order to analyze longer observations we used 
the new 'phase-modulation' or 'sideband' search 
technique described by Ransom et al. (2003) (see 
also Jouteux et al. 2002). This technique relies 
on the fact that if the observation time is longer 
than the orbital period (T b s >> -Porb), the orbit 
phase-modulates the pulsar's spin frequency, re- 
sulting in a family of evenly spaced sidebands in 
the frequency domain, centered around the intrin- 
sic pulsar signal for circular orbits (see for exam- 
ple Figure 2). The width of the region occupied 
by the sidebands is 47r/a;(T bs/^'orb) Fourier bins, 
where / is the frequency of the signal and x is the 
projected semi-major axis. The constant spacing 
between the sidebands is T b s /-P rb Fourier bins. 

A phase modulation search is conducted by tak- 
ing short FFTs of portions of the full power spec- 
trum of an observation. The short FFTs cover 
overlapping portions of various lengths (to account 
for the unknown semi-major axis of the LMXB) 
over all the Fourier frequency range of the origi- 
nal power spectrum (to account for the unknown 
pulsation period). The secondary power spectra 
are then searched for harmonically related peaks 
indicating regularly spaced sidebands (to detect 
the orbital period of the LMXB). Our phase mod- 
ulation search included Fourier interpolation and 
harmonic summing of up to 4 harmonics of the or- 
bital periods. This particular harmonic summing 
accounts for non-circular orbits, but not for non- 
sinusoidal pulsations. A detection provides initial 
estimates of the LMXB spin period, the orbital 
period, and the projected radius of the orbit. For 
interesting candidate oscillations, these estimates 
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are then refined by generating a series of complex- 
valued template responses (i.e. matched filters for 
the Fourier domain response of an orbitally modu- 
lated signal) to correlate with the original Fourier 
amplitudes (Ransom ct al. 2003). 

3.2.2. Obtaining Candidate Pulsation Frequen- 
cies 

To run the phase modulation searches, we first 
concatenated the orbit-long lists of events which 
we had obtained at the end of the preliminary data 
preparation into observation-long lists of events. 
We then binned both the the observation-long and 
the orbit-long lists into time series (see Table 2). 
We performed an FFT on each, and fed the re- 
sults into the local red noise reduction routine. We 
ran the phase modulation searches on the original 
FFTs as well as on the red noise-reduced FFTs. 
The obtained lists of candidate oscillations were 
almost identical: most of the candidates were mil- 
lisecond pulsations so the red noise reduction did 
not affect them. 

In order for the phase modulation search to de- 
tect a phase modulated signal, the short FFTs de- 
scribed in Section 3.2.1 must overlap at least a part 
of the phase modulated signal. When running the 
phase modulation searches on the observation-long 
data sets, we used the range of short FFT lengths 
derived in Ransom ct al. (2003) that was appropri- 
ate for a wide range of orbital parameters: from 64 
to 131072 bins. As we explained in Section 3.2.1, 
the spread of the phase modulated signal in the 
Fourier domain is proportional to xf where x is 
the projected semi-major axis in light-seconds and 
/ is the intrinsic frequency of the signal. Because 
the maximum x for 4U 1820—30 is small (for the 
upper-limit companion mass of 0.3 M Q and an 
inclination angle of 90°, x — 1.02 lt-s), the sizes of 
the short FFTs that we used would allow us to de- 
tect signals up to our Nyquist frequency (2048 Hz). 
We used a smaller range of short FFT sizes when 
we ran the phase-modulation search on the orbit- 
long data sets (64 to 32768 bins), and therefore 
were most sensitive to signal frequencies < 1 kHz 
for all possible values of the companion mass. 

3.2.3. Candidate Follow-up 

For every candidate oscillation it finds, our 
phase modulation search returns, along with other 



information, a pulsation period, an orbital pe- 
riod, and ct. We know from simulation that when 
a phase modulation search detects a signal with 
pulsed fraction equal to or greater than our upper 
limit calculated in Section 4.1, the pulsation pe- 
riod of the signal appears many times in the list of 
candidates obtained for that data set due to over- 
lapping numerous short FFTs. Depending on the 
strength of the signal, there can also be a jump in 
the significance (measured in ct) between the can- 
didate oscillation corresponding to the real signal 
and the candidate with the next smaller signifi- 
cance. When this happens, if the list of candidates 
is arranged in order of decreasing ct, the real signal 
appears at the beginning of the list as an outlier 
with a singularly large ct. 

When examining the phase modulation search 
candidates, we first extracted the candidates 
whose orbital period was 685 s within uncertain- 
ties. Among the extracted candidates, we searched 
for oscillation periods that appeared several times 
in any given observation. There were none. Also 
among the extracted candidates, we searched for 
candidates whose pulsation period appeared in 
the candidate lists of more than 10 observations 
to within a fractional period of 0.1%. There were 
also none. The 0.1% was a generous percentage 
considering that the uncertainties in the candidate 
pulsation periods were much smaller. Finally, we 
searched for outliers within the candidate lists. 
We did not find any. 

3.3. Orbital-Parameter-Fitting Coherent 
Searches 

3.3.1. The Search Technique 

In this search, the source's motion is modelled 
as a circular orbit. This is because LMXB accret- 
ing systems are known to have nearly circular or- 
bits (Zahn 1977). Hence, we ignored the small er- 
rors due to a non-zero eccentricity. Three param- 
eters are required to specify a particular circular 
orbit model: the orbital period P rb, the projected 
semi-major axis in light-seconds x = asin(z)/c 
(where a is the real semi-major axis, i is the in- 
clination of the orbit, and c is the speed of light), 
and the orbital phase at the start of an observa- 
tion = 27r(T sta rt - T )/P 0lh (where T sta rt is the 
start time of the observation, and T a is the time 
of the last ascending node before T sta rt)- Once 
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a set of trial orbital parameters is specified, it is 
straightforward to remove the shifts that the orbit 
would introduce in the arrival time of every pho- 
ton. If the set of trial orbital parameters matches 
the real orbital parameters of the source, the re- 
sulting time series will contain the signature of an 
isolated pulsar that can be found using a standard 
pulsation search with no acceleration. 

For 4U 1820-30, the orbital period is accu- 
rately known, but the other two parameters, x and 
4>, must be searched. In Arons & King (1993), a 
maximum likely value x max for x is estimated, but 
the orbital phase is completely unknown. The sen- 
sitivity of detection to errors in the orbital param- 
eters increases with pulsation frequency (see Sec- 
tion 4.2), so one must also estimate the highest 
likely pulsation frequency / max . From the sensi- 
tivity calculations in Section 4.2 we know that the 
number of sets of trial orbital parameters increases 
as (^max/max) 2 - For 4U 1820—30, reasonable esti- 
mates of these numbers can be provided (see next 
Section), and the effort required to search the data 
is feasible. 

3.3.2. Obtaining Candidate Pulsation Frequen- 
cies 

The method used to run searches by orbital- 
parameter-fitting is essentially the same as the 
standard method for running acceleration searches: 
the arrival times of the photons were time-shifted 
(this time to account for a variable orbital acceler- 
ation), and the resulting data sets were searched 
for isolated pulsar signatures. In particular, 
no matched filtering in the Fourier domain was 
done. The current implementation requires about 
2000 hr to search a typical 4U 1820—30 observa- 
tion on our computers. This is why it was neces- 
sary to limit the number of observations on which 
to run the search. 

Chou & Grindlay (2001) did a folded period 
searching of the combined 1996-1997 RXTE stan- 
dard2 observations of 4U 1820—30. They obtained 
an orbital period of P orb = 685.0144 ± 0.0054 s. 
They then did a phase analysis combining the 
RXTE 1996-1997 observations of 4U 1820-30 
with pre-XTE data and found an orbital period 
of P O rb=685.0119 ± 0.0001 s. The uncertainties 
in both figures are small enough (see Section 4.2) 
that we did not need to search over this parame- 
ter. The decay in P or b due to emission of gravita- 



tional radiation, estimated to be ~ 0.00018 s yr _1 
(Wagoner 1975), fell within the above uncertain- 
ties. We used P or b = 685.014 s when we ran our 
orbital-paramctcr-fitting search on 4U 1820—30. 

The maximum plausible x, assuming the com- 
panion is a degenerate star is 0.024 lt-s (this 
corresponds to a companion mass of 0.078M Q 
for an edge-on inclination). With our orbital- 
parameter-fitting search, we searched up to an 
i ma i of 0.032 lt-s (corresponding to a compan- 
ion mass of 0.106M Q for an edge-on inclination. 
Larger companion masses are allowed for smaller 
inclinations). We also set an upper bound on 
searched pulsation frequency of / max = 1250 Hz 
to save computer time. If a set of trial orbital 
parameters corresponds exactly to the real orbital 
parameters of 4U 1820—30, our orbital-parameter- 
fitting search collects successfully all the power 
due to the signal. If our trial orbital parame- 
ters are slightly off, there will be a sensitivity loss. 
When we ran the search, we tuned the parameter 
steps to accept up to 50% loss in signal power at 
the maximum frequency. (As frequency decreases 
or the error in the size of the orbit decreases, this 
loss decreases to very near zero; see Section 4.2). 

Based on calculations by in Middleditch (2004) , 
and using the above known properties of 4U 1820—30 
and the above mentioned maximum loss in signal 
power, we calculated that the error in a trial x 
should be no more than 0.000143 lt-s (see Section 
4.2 for details). With that in mind, we generated 
an optimal grid of trial orbital parameter sets. 
The grid contained several thousands of points 
(see Section 4.2) 

For each set of trial orbital parameters, we ran a 
separate search. We time-shifted the event data, 
to reflect the trial orbital parameters. We then 
binned it at a resolution of 0.244 ms and com- 
puted its Fourier transform. We applied our lo- 
cal red-noise reduction code and fed the result to 
our standard acceleration search code with a max- 
imum acceleration of zero and an increased sig- 
nificance threshold (i.e. a > 3); keeping in mind 
that the a returned by the acceleration search af- 
ter subtracting the time-shifts is an overestimate 
of the real a of the signal. This is because the 
orbital-parametcr-fitting coherent search fails to 
take into account the number of times accelera- 
tions searches were run. The resulting list of can- 
didate oscillations was then stored and the process 
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was repeated for the next set of trial orbital pa- 
rameters. 

3.3.3. Candidate Follow-up 

Each candidate returned by the orbital-parameter- 
fitting code included a pulsation frequency /, an 
orbital amplitude x, an orbital phase (f> 7 and a 
significance a. Using the returned values for /, x, 
and (j), we folded all the candidates with a > 5, 
allowing the folding program to search a small dis- 
tance in the frequency space. Most of the folded 
candidates showed noise-like profiles. None of the 
candidates corresponding to the more promising- 
looking folds appeared in more that one obser- 
vation except for one candidate. This candidate 
appeared in two observations. However, the fold 
of the second detection of this candidate showed 
a noise-like profile. 

Since folding the statistically significant candi- 
dates did not return any interesting results, we 
looked for candidates that repeated in more than 
one observation (to within a fractional range of 
/percent = 0.001% and to within a difference in x of 

•^•difference — (-Porb * /percent)/ (4007T + 27r/ percen t) lt- 
s)« /percent 

^difference were derived taking into 
account the spacing between the trial values of x 
and the drop in the amplitude of the detected sig- 
nal when the trial orbital parameters were slightly 
off. There were about 300 candidates that re- 
peated in 2 observations, and candidates that re- 
peated in more than 2 observations. We grouped 
the 300 candidates into pairs. Assuming the or- 
bital period of 685.0119 ± 0.0001 s reported in 
Chou & Grindlay (2001), we checked if the orbits 
of the members of each pair were in phase with 
each other. In 18 of the pairs the two orbits were 
in phase. In all 18 cases the single trial signifi- 
cances of the members of the pairs were around 
3cr. Again, the folds of all 36 candidates showed 
noise-like profiles. 

4. Results and Upper Limits 

Our searches did not detect pulsations. 

Since we ran the acceleration searches in a re- 
gion where they are not optimized, they were 
much less sensitive than the other two searches 
that we used. Below we explain the procedure 
we used to set an upper limit on the pulsed frac- 
tion of pulsations from 4U 1820—30 detectable 



with our phase modulation search technique and 
our orbital-parameter-fitting coherent search tech- 
nique. 

4.1. Phase Modulation Searches 

The phase modulation searches that we ran 
on observation-long data sets were more sensitive 
than the searches that we ran on the orbit-long 
data sets (Ransom et al. 2003). In this section we 
present the procedure used to set an upper limit on 
the pulsed fraction using the more sensitive search. 

First, to save computer time, we wrote a modi- 
fied version of our search script that searched only 
a particular region of the FFT. We then ran this 
version of the search on several data sets con- 
taining simulated pulsations from 4U 1820—30 of 
various pulsed fractions, for a range of compan- 
ion masses (0.02 M Q < M c < 0.3 M Q ), and for 
several spin frequencies. For the purpose of the 
simulations, we assumed the known orbital pe- 
riod of 685 s, a typical RXTE observation duration 
(20.5 ks), a count rate of 2050 counts per second, 
an inclination angle of 60°, and a circular orbit. 
Figure 3 shows the results of the simulations for 
this type of search run on data sets containing fake 
3 ms sinusoidal pulsations. Every point in the fig- 
ure corresponds to 100 searched data sets. The 
shading of every point corresponds to the detected 
power of the signal. The power was then converted 
to probability and multiplied by the number of in- 
dependent trials (roughly the original number of 
time bins) to convert it to a. In Figure 3, the line 
of a = 5 indicates the value of our 95% confidence 
upper limit on the pulsed fraction for a signal pe- 
riod of 3 ms. This value varied between 0.55% and 
1.0% (background unsubtracted) over the range 
of companion masses that we used. The value 
was unaffected by the background correction since 
the background was small compared to the source 
count rate: using the 'FTOOL' PCABACKEST, we 
estimated the background count to be 10 counts 
per second per PCU. This is only 2% of the typ- 
ical count rate of the source. This means that 
if there were a 3-ms coherent pulsation coming 
from 4U 1820—30 with a pulsed fraction equal to 
or higher than the upper limit stated above, our 
searches would have detected it at a > 5 95% of 
the time. The upper limit on the pulsed fraction 
decreased slightly when the period of the signal 
was larger by a few milliseconds. This is because, 
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at longer spin periods, the width of the sidebands 
feature is smaller (see Section 3.2.1); and hence 
the power accumulated in each of the sidebands is 
greater. For example, for a 10-ms coherent pul- 
sation, this value varied between 0.3% and 0.5% 
over the range of companion masses that we used. 
Similarly, the upper limit on the pulsed fraction 
increases for smaller spin frequencies. 

Our upper limit on the pulsed fraction of pos- 
sible signals from 4U 1820-30 is -0.8% for spin 
periods close to 3 ms. The pulsed fractions of the 
signals detected from the 5 LMXBs with known 
signal period ranged between 3% and 7%, for sig- 
nal frequencies between 2 and 6 ms (Wijnands & 
van der Klis 1998; Markwardt et al. 2002; Gal- 
loway et al. 2002; Markwardt et al. 2003a; Cam- 
pana et al. 2003). 

4.2. Orbital-Parameter-Fitting Coherent 
Searches 

The goal of orbital-parameter-fitting is to 
achieve a sensitivity close to the sensitivity one 
would have if the pulsar were isolated. To achieve 
this, one must compute how densely to sample 
the three-dimensional space of orbital parame- 
ters. Such calculations have already been done by 
Middleditch (2004), but in this section we work 
through them with more details to better under- 
stand the effects of small errors in the trial orbital 
period. 

First, assume that P or b is known exactly. Then 
combine x and <j) into a single complex number X 

X = xe*+. 

If the pulsar emits sinusoidal pulsations at a fre- 
quency /, then the received signal will be 

s{t) = sin(27r/ (t + Re (x orb e 2 ™'/ p " b ))), 

up to a phase shift (which we will ignore). If the 
pulse profile is not sinusoidal, this transformation 
will happen equally to all harmonics; we will ig- 
nore this too. 

Suppose that we remove the effects of the or- 
bital motion by time-shifting the incoming events 
according to a model using the known P or b but 
X = -X" tr iai- The observed signal will then be 

s(t) = sin(27r/(i + 



Re (x orb e 27 " t/P " b ) - Re (x tri aie 27H ' /P " b ))) 
= sin(27r/(i + 

i?e((X orb -X trial )e 2 ^/ p » b ))). 

In other words, the effect will be exactly that of a 
residual phase modulation of complex amplitude 

X rb — -Atrial ■ 

Such a phase modulation produces a comb- 
like pattern in the Fourier domain (see Figure 2). 
For narrow combs, the tallest peak is the cen- 
tral peak: the one at the true pulsation fre- 
quency The amplitude of the central peak falls 
off as JoG-^orb — -Atrial | 27t/) where J is a th 
order Bessel function (Ransom et al. 2003; Mid- 
dleditch 2004). Therefore, the sensitivity loss (in 
power) due to a residual phase modulation of size 
\X or b — trial I = £ will be Jo(s2tt f) 2 . This means 
that if we search the X OI b space with a sufficient 
density to fall within e of every plausible X OI ^, 
we will detect the source with a sensitivity only a 
factor of J (e27r/) 2 less than that we would have 
for an isolated pulsar. In other words, the resid- 
ual time shifts need to be small compared to the 
pulsation period. 

For this particular situation, we fixed the maxi- 
mum \X\ as 0.032 lt-ms, / max as 1250 Hz, and the 
sensitivity factor as 0.5 at f max . The maximum 
allowable residual X (which we called e) is then 
0.000143 ms. With that in mind, we generated 
an optimal grid of trials for complex- valued X. It 
turns out that the optimal grid (the one where the 
parameter space is searched most efficiently) is a 
a circle-packing array: a grid where the values of 
-Atrial are distributed on the corners of equilateral 
triangles (e.g., Middleditch 2004). The point in 
the middle of every triangle is the point that is fur- 
thest away from the closest X tr i a i and therefore the 
distance between it and the closest -X" tr i a i should 
be e — 0.000143 lt-s. The spacing between any 
two guesses is then 0.000264 lt-s. The total num- 
ber of trials of X is then the area of the complex 
X-plane containing the trials (roughly 7ra; max 2 ) di- 
vided by the area covered by every X tr i a i (roughly 
7T£ 2 ). This amounts to an approximate number 
of 50000 trials per observation. With the trials 
of X distributed as we described, our worst case 
sensitivity to a given signal of frequency / is then 
Jo(0. 000143 x 27r/) 2 x (the sensitivity to an iso- 
lated pulsar signal of frequency /). At / = / max 
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the value of this factor is 0.5. Its value increases 
for signals of lower frequency. For / = / ma x/2, it 
is 0.85. 

In the above calculation, we assumed no loss of 
power due to errors in the assumed orbital period. 
But how accurately does P or b need to be known? 
We can do an approximate calculation: suppose 
that Ptriai = (1 — e)-P rb- The resulting signal is 
then 

s(t) = sin(27r.f(t + 

Re (x orb e 2 -*/ p -) - Re (x triaie 2 ^-) )) 
~ sin(27r/(t + 

Re ((A orb - X tria ie 2 " t£/p » b ) e 2mt/p - b ))). 

This amounts to an additional slowly-varying er- 
ror in the estimation of the orbital phase <f). Errors 
in that are independent of (P or b — -Ptriai) have al- 
ready been taken into account by the earlier sensi- 
tivity calculations. Therefore, when we calculated 
the effects of the additional slowly- varying error in 
4> in a given observation of length T b s , we could 
assume, to start with, that (j> is correct in the mid- 
dle of that observation. At the end of that obser- 
vation, we would have accumulated an error of size 
e7rT bs/-Porb in the estimation of 4>. This adds an 
additional residual X of about e \X\ 7rT b s /Porb- 
In the worst case, this is about 4.3e lt-s. From the 
difference between the value we used for the orbital 
period (685.014 s) and the lower bound on the 
more accurate orbital period measured in Chou & 
Grindlay (2001) of 685.01 19±0. 0001 s, we estimate 
e ~ 0.000003. Hence the size of the maximum ad- 
ditional residual X is ~ 0.000013 lt-s, which is a 
factor of ~11 smaller than the maximum residual 
X that the spacing in our grid of trials is allow- 
ing. In the worst case scenario (at / = / max and 
x = :Cmax), if we combine the maximum error in X 
due to the spacing of our trials with the additional 
error in X due to the error in the orbital period, 
we find a worst case sensitivity of 0.43 x that for 
an isolated pulsar (instead of 0.5 when the P or b is 
exactly right). If the error caused by the uncer- 
tainty in P or b were slightly larger, it would have 
caused an even larger decrease in the detectable 
power, and we would have been forced to search 
in P orb also. 

Above, we have extracted formulae for the sen- 
sitivity of the orbital-parameter-fitting search in 



terms of the sensitivity of searches for an isolated 
pulsar. These, as well as a small number of simula- 
tions, show the following: if there are no errors in 
the search parameters, we can detect at 3 a a sig- 
nal of pulsed fraction 0.17%, and at 5.3 cr a signal 
of pulsed fraction 0.2%. If our trials for the orbital 
parameters (including the orbital period) are off 
by the worst-case amounts, we can detect at 3 a a 
signal of pulsed fraction 0.26%, and at 5.3 cr a sig- 
nal of pulsed fraction 0.3%. The above numbers 
are summarized in Table 3. Because of the way 
we are doing the candidate follow-up (i.e. folding 
candidates that appear several times in the list of 
candidates), the occurrence of the signal in more 
than one observation increases greatly the chance 
of detecting it: if a signal appeared only in one 
observation at 3 cr, we would have discarded it as 
a false candidate. If a signal appeared in several 
observations at 3 a we would certainly not have 
missed it. Finally, any individual signal that ap- 
peared at a > 5.3 would have been detected with 
confidence. 

To summarize, our sensitivity formulae and our 
simulations show reliable detection at a pulse frac- 
tion of 0.3% for the parameters we have chosen. 
This is of course dependent on a number of as- 
sumptions: the companion must be a degenerate 
star, the pulsation frequency must be less than 
about 1250 Hz, and P or b must be within 0.0022 ms 
of 685.014 ms. It should be noted, however, that 
these assumptions are covariant: if the pulsation 
frequency is lower, the sensitivity improves. If x 
is smaller than the assumed upper bound, we be- 
come less sensitive to an error in Ptriai- 

5. Discussion 

The search described in this paper represents 
the first source in a larger program to search for co- 
herent pulsations in the persistent LMXBs. There 
are currently five known accreting MSPs, all of 
which have pulse frequencies in the range of 185 
to 435 Hz and low luminosities. The pulsed frac- 
tions observed are in the 3 - 7 % range. Recently, 
observations of burst oscillations, kHz QPOs, and 
coherent pulsations in two of these accreting MSPs 
have confirmed that the burst oscillations are in- 
deed at the spin frequency (Chakrabarty et al. 
2003; Markwardt, Strohmayer, & Swank 2003b) 
and that the kHz QPO separation can be ei- 
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ther near the spin frequency or half the spin fre- 
quency. Thus, there is now a total sample of 
14 LMXBs with well-determined spin frequen- 
cies. In 5 LMXBs, coherent pulsations are ob- 
served (Wijnands & van der Klis 1998; Markwardt 
ct al. 2002; Galloway et al. 2002; Markwardt et al. 
2003a; Markwardt & Swank 2003). In 11 LMXBs 
burst oscillations are observed (Strohmayer 2001; 
Strohmayer et al. 2003). In 2 of these LMXBs, 
both coherent pulsations and burst oscillations are 
observed (Chakrabarty et al. 2003; Strohmayer 
et al. 2003). 

We can see pulsations during X-ray bursts and 
in the non-burst emission in a handful of transient 
sources. Why then can we not see pulsations in 
the persistent emission of the many bright LMXBs 
including 4U 1820-30? 

For any individual pulsar, an unfavorable view- 
ing geometry or rotational geometry can prevent 
the pulsations from being seen. However, a sys- 
tematic absence of pulsations from all persistent 
LMXBs would require a stronger explanation. 
There are three main categories of explanations 
that have been considered. These are: relativistic 
effects that decrease the observed pulsed fraction, 
smearing of pulses from the surface due to elec- 
tron scattering, and surface magnetic fields that 
are too weak to significantly channel the accretion 
flow (possibly as a result of magnetic screening). 
We briefly consider each of the possibilities below. 

In Wood, Ftaclas, & Kearney (1988), the 
authors demonstrate that gravitational light- 
bending effects on emission from hot polar caps 
can greatly reduce the pulsed fraction: the pulses 
are distorted and the pulsed fraction decreases as 
a star becomes more relativistic because its radius 
is decreasing. The decrease in the pulsed frac- 
tion results from the fact that the gravitational 
deflection of photons displaces the horizon away 
from the observer. In other words, more of the 
star than a hemisphere is visible to the observer. 
Thus, one of the caps is always in nearly full view 
as the other is disappearing. It is possible for the 
pulsed fraction of a given source to be reduced 
down to ~1%, provided the radius of the emit- 
ting star is small enough and the opening angle 
of the polar caps is large enough (Meszaros, Rif- 
fert, & Berthiaume 1988). This fraction is larger 
than our upper limit on the pulsed fraction, sug- 
gesting that gravitational lensing effects may not 



be the only factor preventing us from seeing the 
pulsations. Also, the five known accreting MSPs 
exhibit pulsed fractions of more than 5%. There 
is no reason to believe that those systems have a 
significantly different value of M/R than the per- 
sistent LMXBs. Consequently, it seems that the 
relativistic effects are not the primary explanation 
for the small pulsed fraction in the latter sources. 

Another explanation for low modulation frac- 
tions, electron scattering, has been revisited by 
Titarchuk, Cui, & Wood (2002). The existence 
of a corona (a spherical cloud) of hot electrons 
surrounding the accreting star has been inferred 
from the analysis of the X-ray spectra of several 
LMXBs (e.g., Kylafis & Klimis 1987; Titarchuk 
ct al. 2002). This hot corona is presumably up- 
scattering softer seed photons and smearing out 
the pulsations. Titarchuk et al. (2002) argue that 
it is because of the large optical depth of this 
corona that pulsations have not been detected 
from most LMXBs. They claim that the opti- 
cal depth below which millisecond pulsations be- 
come detectable is r ~ 4. They also do spectral 
fits for several LMXBs and extract the parame- 
ter r. For the majority of the analysed LMXBs 
they find that the optical depth of the Compton 
cloud is at least twice as high as t = 4. For these 
LMXBs, they claim that any high frequency pul- 
sations (higher than ^300 Hz) would be below de- 
tectability. Titarchuk et al. also make the obser- 
vation that there is a possibility of finding pulsa- 
tions in the sources that are much less luminous, 
because in these sources, the spherical Compton 
cloud is more transparent (r < 4). Unfortunately, 
the smallest reported r for the spherical Compton 
cloud that is presumably surrounding 4U 1820—30 
was t = 6 (Bloser et al. 2000). In all the other 
papers that discuss spectra from 4U 1820—30, the 
authors used a Comptonization model with a black 
body component to fit the spectra and found that 
r falls between 11 and 13 (e.g., White, Stella, & 
Parmar 1988). If what Titarchuk et al. (2002) 
have claimed is true, the large optical depth of the 
Comptonizing cloud could very well be the reason 
we are not seeing pulsations from 4U 1820—30. 

The possibility of magnetic screening by ac- 
creted material at high accretion rates has been 
examined by Cumming, Zweibel, & Bildsten 
(2001). Their reasoning goes as follows: to get 
preferential heating of the stellar surface, the mag- 
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netic field must be large enough to disrupt the 
accretion flow and channel it onto the polar caps. 
Preferential heating of the surface causes the ob- 
served emission to be pulsed. If the accretion flow 
is not disrupted by a magnetosphere, the accre- 
tion disk is believed to extend almost to the stellar 
surface and pulsations arc no longer seen. This 
could happen because the external dipolc mag- 
netic field of the star is 'buried' or 'screened' by 
the accreted matter. Cumming et al. (2001) study 
the 'burial' scenario. In particular, they study 
the competition between downward compression 
of the magnetic flux by accretion at a given M, 
and its upward transport by diffusion. They find 
that the diffusion and accretion time scales are 
equal for M ~ O.lMEdd, where MEdd is the ac- 
cretion rate given by the Eddington limit. This 
means that for larger M there is little chance that 
magnetic flux makes it to the surface and there- 
fore little chance of observing pulsations. They 
also find that magnetic screening is completely 
ineffective for M < 0.02MEdd, so that, no mat- 
ter how the accreted material joins onto the star, 
the underlying stellar field would always be evi- 
dent. The typical range of mass accretion rates in 
LMXBs is M > 0.01M E dd (M > 10~ 10 M yr" 1 ) 
(Chakrabarty et al. 2003). Cumming et al. (2001) 
propose that the absence of persistent millisecond 
pulsations from most neutron stars in LMXBs 
is due to their high accretion rates. Discoveries 
made since then are consistent with their hypoth- 
esis: all five persistent millisecond X-ray pulsars 
are low-luminosity transients in very close binaries 
that lie at the low end of the mean M distribution 
for LMXBs with M ~ 10" 11 M yr" 1 (Galloway 
et al. 2002; Campana et al. 2003; Chakrabarty 
et al. 2003). It is tempting to adopt magnetic 
screening as the main explanation for the lack of 
observed pulsations from 4U 1820—30 because the 
luminosity of this source (~ 4 x 10 37 ergs/s) is on 
the high end of luminosities of atoll sources, imply- 
ing a normal (rather than a low) accretion rate of 
the in-falling matter. Stella et al. (1987) estimated 
this accretion rate to be ~ 5 x 10~ 9 M Q yr" 1 . 

6. Summary 

In this paper, we have described a search of 
4U1820— 30 that set upper limits on the pulsed 
fraction covering essentially all plausible pulse 
periods and companion masses for this system. 



These limits were between 0.3% and 0.8% over this 
whole range. We also searched a broad range of 
luminosity and spectral states (see Figure 1). The 
luminosities covered the range from 0.08 to 0.25 
of the Eddington luminosity for a 1.4 M Q neu- 
tron star and a distance of 6.4 kpc. 4U1820— 30 is 
interesting in this regard because the luminosity 
varies by a factor of three. 

For 4U 1820—30, our upper limits are about an 
order of magnitude lower in pulsed fraction than 
the known accreting MSPs. However, this is only 
one source and any number of specific features, 
such as unfavorable viewing or rotational geome- 
tries, in addition to electron scattering and mag- 
netic screening, could prevent the pulsations from 
being seen. The most interesting conclusions will 
be drawn once we have a sample of many more 
sources that have been searched to similarly deep 
levels. 
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Fig. 1. — The power-law spectral index and the 2- 
10 keV luminosity of each of the 34 observations 
that we have searched for pulsations. These num- 
bers were obtained from spectral fittings and were 
based on an assumed distance of 6.4 kpc. 
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Period - 0.00229717150807573 (s) 
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Frequency - 435.317955047619 (Hz) 

Fig. 2. — Sidebands in the FFT of observa- 
tion 70131-01-05-00 (4 RXTE orbits long) of 
XTE J 175 1-305 centered on the intrinsic pulsar 
spin frequency (435.318 Hz) and caused by or- 
bital phase modulations of the pulsar frequency 
(see Section 3.2.1). The binary orbital period that 
is causing the above modulation is ~42.4 min and 
the projected orbital radius is ^10.11 lt-ms. The 
pulsed fraction above is high (3.76%, background 
unsubtracted) which makes the signal stand well 
above the noise. The width of the sideband pat- 
tern is related to the projected semi-major axis of 
the orbit. The above signal was detected at 10.6 a 
using our phase modulation search technique and 
at 10.8cr using our acceleration search technique. 
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Fig. 3. — Monte-Carlo derived sensitivity calcula- 
tions for the pulsation searches of archival RXTE 
observations of 4U 1820—30. The plotted sensi- 
tivities are 95% confidence limits using the phase 
modulation search technique for a typical RXTE 
observation assuming a pulsar spin period of 3 ms. 
With these types of searches, we would easily have 
detected any coherent pulsations with a pulsed 
fraction > 0.8% for all realistic companion masses 
and spin periods. 
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Table 1 

RXTE Observations' 1 of 4U 1820-30 





ObsID 




Year 


Start Timc b 


Duration c 


Time on 


Count Ratc d 


Num 


Num Events f 


Luminosi 










(M.TD) 


(ks) 


Source (ks) 


(counts/s/PCU) 


PCUs c 


(X 10 6 ) 


(X 10 37 erg 


1 


10075-01-01- 


020 


1996 


50384.63103 


28.0 


18.2 


450 


5 


38.56 


2.4 


2* 


10075-01-01- 


02 


1996 


50384.96433 


17.2 


10.4 


460 


5 


26.08 


2.5 


3 


10075-01-01- 


■030 


1996 


50386.49418 


28.2 


17.1 


600 - 565 


5 


44.99 


3.0 


4* 


10075-01-01- 


031 


1996 


50386.82749 


25.2 


16.4 


565 - 455 


4 


42.17 


2.7 


5 


10075-01-01- 


03 


1996 


50387.14749 


3.1 


3.1 


480 


4 


7.50 


2.6 


6 


20075-01-01- 


00 


1997 


50488.42272 


24.9 


16.0 


635 - 550 


5 


48.30 


3.3 


7 


20075-01-02- 


■01 


1997 


50514.02219 


17.9 


11.1 


725 


5 


42.19 


3.8 


8 


20075-01-03- 


00 


1997 


50530.70491 


20.0 


13.1 


585 


3 


33.96 


3.0 


9* 


20075-01-04- 


00 


1997 


50548.70709 


24.1 


14.9 


450 


5 


33.98 


2.4 


10* 


20075-01-06- 


■00 


1997 


50595.33191 


15.3 


11.3 


720 


5 


41.23 


3.9 


11* 


20075-01-08- 


00 


1997 


50645.22308 


15.2 


11.0 


720 


5 


40.94 


3.8 


12* 


20075-01-09- 


■00 


1997 


50675.45754 


25.1 


16.0 


670 - 735 


4 


51.31 


3.6 


13 


30057-01-01- 


■ 02 


1998 


50909.66291 


18.4 


11.2 


470 


4 


24.24 


2.5 


14* 


30053-03-01- 


■000 


1998 


50920.19895 


25.2 


16.4 


465 




39.94 


2.5 


15 


30053-03-01- 


■001 


1998 


50920.52779 


20.8 


12.0 


465 


4 


24.94 


2.5 


16* 


30053-03-01- 


■00 


1998 


50920.79458 


25.8 


13.2 


475 




41.05 


2.6 


17 


30053-03-02- 


■03 


1998 


51080.45284 


20.9 


14.2 


430 - 380 


4 


22.30 


2.2 


18 


40017-01-02- 


■00 


1999 


51222.63887 


20.9 


14.4 


510 


4 


33.63 


2.9 


19 


40017-01-03- 


00 


1999 


51238.03748 


15.1 


10.4 


320 


5 


16.81 


1.8 


20 


40017-01-04- 


■00 


1999 


51253.69624 


18.9 


10.5 


430 


4 


13.42 


2.4 


21* 


40017-01-07- 


■00 


1999 


51300.07107 


18.6 


10.9 


980 - 895 


3 


37.42 


5.0 


22 


40017-01-09- 


■00 


1999 


51330.51073 


15.1 


10.6 


700 


3 


28.80 


3.9 


23 


40017-01-12- 


■00 


1999 


51389.38952 


15.3 


10.7 


705 


3 


23.53 


3.8 


24* 


40017-01-13- 


00 


1999 


51400.44912 


21.0 


14.7 


660 


3 


28.88 


3.6 


25 


40017-01-14- 


■00 


1999 


51407.50805 


15.1 


10.6 


505 


5 


27.13 


2.9 


26 


30057-01-03- 


■01 


1999 


51410.50430 


24.0 


15.1 


515 - 445 


3 


24.59 


2.5 


27 


30057-01-03- 


-000 


1999 


51412.31866 


29.0 


16.9 


495 - 530 


3 


34.10 


2.9 


28 


30057-01-03- 


■00 


1999 


51412.65197 


7.5 


5.1 


510 


4 


10.74 


2.7 


29 


40017-01-15- 


■00 


1999 


51417.29843 


20.4 


12.8 


440 


3 


17.05 


2.4 


30 


40019-02-01- 


■00 


1999 


51421.30667 


25.0 


15.3 


370 


3 


16.81 


2.0 


31 


30057-01-05- 


00 


1999 


51435.00244 


26.8 


17.9 


365 - 480 


2 


26.34 


2.4 


32 


40017-01-18- 


■00 


1999 


51464.36730 


20.8 


14.0 


535 


3 


26.07 


3.0 


33 


40017-01-21- 


■00 


2001 


51941.01232 


14.9 


10.4 


380 


4 


14.36 


2.0 


34* 


40017-01-23- 


■01 


2002 


52355.88359 


14.8 


10.3 


915 


3 


32.89 


4.8 



Wc used event data with ~125 /J.s time resolution for all observations. 
b Start time of the observation in MJD. 
c Stop time — Start time of the observation in ks. 

d Approximate photon count per second per PCU obtained from the Standard 2 lightcurve in the wide band (2-20 keV). 
°Minimum number of PCUs in use at a given moment of the observation. 
^Approximate total number of events processed in the wide band (2-20 keV). 

s 2-10 keV luminosities obtained from spectral fitting and based on an assumed distance of 6.4 kpc. 

* Observations marked with a * were searched with acceleration searches, phase modulation searches, and an orbital-parameter-fitting coherent search. 
The remaining observations were searched only with acceleration searches and phase modulation searches. 
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Table 2 

Search Types for Different Data File Types 



Type of Data File 



Searches Run 



Observation-long lists of events 

Observation-long time scries binned at 0.244 ms 
KXT&orbit-long lists of events 

fiXTE-orbit-long time series binned at 0.122 ms 
600-s-long and 1500-s-long lists of events 



orbital-parameter- fitting coherent search 3, (W band only) 

Acceleration searches (all bands), phase modulation searches (all bands) 
No searches run 

Acceleration searches (all bands), phase modulation searches (all bands) 
No searches run 



600-s-long and 1500-s-long time series binned at 0.122 ms Acceleration searches {all bands) 



a Bccausc running an orbital-parameter- fitting coherent search was extremely demanding computationally, we ran this search on event lists of 11 
observations only. By contrast, we ran accelerations searches and phase modulation searches on time series of all 34 observations. 

b Thc segmented acceleration searches that we ran, even though they were not short enough for the acceleration to be constant, could still be 
useful if a strong pulsed signal were to appear for a short period of time. 
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Table 3 

Results of Sensitivity Calculations for the 
Orbital-Parameter-Fitting Coherent Searches 



Pulsed Fraction 


Best Case 


Worst Case 




Detection a* 


Detection a* 


0.17% 


3 


< 3 


0.2% 


5.3 


< 3 


0.26% 


> 5.3 


3 


0.3% 


> 5.3 


5.3 



*er is an equivalent Gaussian significance. 
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